Two-fermion bound state in a Bose-Einstein condensate 
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A nonlinear Schrodinger equation is derived for the dynamics of a beam of ultracold fermionic 
atoms traversing a Bose-Einstein condensate. The condensate phonon modes are shown to provide 
a nonlinear medium for the fermionic atoms. A two-fermion bound state is predicted to arise, and 
the signature of the bound state in a nonlinear atom optics experiment is discussed. 
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The recent experimental success in creating quantum- 
degenerate atomic Fermi gases [0, ||, ||, [|, [| is opening 
up fascinating new opportunities to explore the quantum 
statistics of ultracold atoms. Of particular interest in this 
context is the formation of atomic Cooper pairs in ultra- 
cold Fermi gases ff, g, §, §, Q 0, |l|. Another impor- 
tant question is whether nonlinear atom optical effects, 
which are now well established in bosonic systems, can 
also occur with fermions. Fermionic behavior is strongly 
constrained by the Pauli exclusion principle. This limits 
the variety of possible nonlinear atom optics effects, but 
also offers the potential for novel applications without 
analogs in optics. These include, for example, low-noise 
inertial and rotation sensors, and quantum information 
processing. Research on nonlinear atom optical effects in 
a quantum-degenerate Fermi gas has recently been ini- 
tiated. Two theoretical papers |l4| have indicated 
that fermionic atom four-wave mixing is possible under 
appropriate conditions. 

Our goal in this Letter is to explore a new situation 
where a gas of bosons serves as a nonlinear medium for 
fermionic atoms. In particular, we study how the inter- 
atomic interaction between a Bose-Einstein condensate 
and a fermionic beam can be employed to manipulate the 
quantum state of the beam. By drawing on the analogy 
to nonlinear optics we describe the interaction in terms of 
an effective attractive Kerr nonlinearity, and show that 
a two-fermion bound state can result with a unique sig- 
nature in a nonlinear atom optical experiment. 

Since the Pauli exclusion principle precludes the direct 
evaporative cooling of spin-polarized fermionic samples, 
current experiments employ either unpolarized fermionic 
mixtures |IJ U or Bosc- Fermi mixtures ||, ||, |[. Our 
starting point is a beam of fermionic atoms with two 
internal spin states | f) and | J.) interacting via two- 
body collisions with a quantum degenerate Bose gas. The 
atomic Bose field is decomposed in the familiar way as 
4> B (r) = </>s(r) + £(r) = \fn B {r) + £(r), where B (r) is 
the condensate wave function, taken to be real for sim- 
plicity, £(r) describes the elementary excitations above 
the condensate, and tib = \4>b\ 2 is the condensate den- 
sity. The fermionic field is described by the field oper- 
ators V'cr(r), with a = {t,l}- In the following, it will 
also be useful to introduce the density fluctuations of the 



fermion fields by 5h a — ?/>£(r)?/v(r) — (r)^ CT (r)) = 
h<j - n a . 

The Bose and Fermi systems are coupled by two-body 
interactions. In the shapeless approximation, and taking 
into account that s-wave scattering is forbidden between 
fermionic atoms of same spin, we find that to lowest- 
order in £(r), the dynamics of the bosonic atoms is given 
by the coupled equations Jll| 

To the same order, the fermionic field equations are 

+ +4> b £)^+w14<t4*- (2) 

Here, H { ° ] =f a + V a (a = B,a) are the single-particle 
Hamiltonians for bosonic atoms and for fermionic atoms 
of spin cr, respectively, T a and V a are the associated 
kinetic energy and trapping potential. The Hamilto- 
nians H B = + 2gn B + $e n o> an d H Fa = 
Hp a + fcrUB, also include the self-contribution to the 
mean-field energy of the respective fields. Finally, fi a 
are the chemical potentials, and the parameters g, f 
and h represent the boson-boson, boson-fermion, and 
fcrmion-fcrmion interaction strengths. In terms of the 
s-wave scattering lengths a, they are given by g — 
Airh 2 (lb /rriB, fa — 2nh 2 aBF<7 /m r , h = ATifi 2 a^i/mF, 
with m r — m pm b / {ttlf + tub) being the reduced mass. 

In previous work JTl[ we discussed how the coupling to 
a fermionic component can change the dynamical stabil- 
ity of a Bose condensate. Here we concentrate instead on 
how the presence of a condensate can induce nonlinear 
dynamics of a fermion field. To isolate the key underlying 
physical mechanisms, it is useful to simplify the situation 
as much as possible. With this in mind, we consider a 
situation where the back-action of the fermionic fields on 
the condensate is negligible. Specifically, we assume that 
griB >> J2a f<? n °- ls satisfied, in which case we can ig- 
nore the effects of the fermionic beam on the condensate 



2 



wave function 0b (r). We then apply a standard Bogoli- on the excitation field £(r). One readily finds 
ubov approach [jl5| to determine the effect of the fermions 



Kr, t) = C (0) (r, t) + 1 jT dt' | dV [G(r, r', t - 0<MO - F(r, r', i - 0^(0] E f' 6 **- 



(3) 
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The first term ^°- ) (r, t) on the right-hand side describes 
the free-field vacuum quasi-particle fluctuations in the 
absence of fermions. It has the familiar form 



|(°)(r,i) 



(un(r)e 1 



-iE n t/H, 



<(r)e 



>,E n t/h 



with the Bogoliubov quasi-particle operators a n and 
ajj satisfying Bose commutation relations. The quasi- 
particle mode functions u n {r) and v n (r), and correspond- 
ing energy eigenvalues are determined by the matrix 
equations 



gn B 



-gn B 
-Hb 



Mn(r) 



E n 



u n (r) 



(4) 



The second term in Eq. (^|) is a four-wave mixing pro- 



cess that mixes the condensate with quasi-particles. It is 
mediated by the density fluctuations of the Fermi fields, 
whose evolution is governed by the quasi-particle Green's 
functions 

G(r, r', r) = E (e^^u^u^r') - e^<(r)i; n (r')J 



with a similar form for F (r, r',r), but with u'(r') and 
v n (r') replaced by (r') and u a (r'), respectively. 

The lowest-order contribution of the condensate to the 
dynamics of the Fermi fields is obtained by substituting 
Eq. (||) into Eq. (||) . Using the homogeneous case quasi- 
particle mode functions of Ref. for simplicity, this 
yields the stochastic Heisenberg equations of motion 



$ a + Wl,i><x>i>o+ V] f d 3 r' [ dTW aal (r,r',T)(j) B (r)^ B {r')dh a ,{r',t-T)ip„ + f (T ip a . (5) 



r 



Here we have defined 



W aa , (r, r', r) = ( - ) [A m , (r, r', r) - A; ff , (r, r', r)] , 



fcrfc 



V 



3 -iB k T/fi+ik-(r— r') 



with ek = h 2 k 2 /2rriB, and -Ek = ^/ £k(ek + 2gns). 

The third term in Eq. (^|) is the nonlinear optics ana- 
logue of a fifth-order nonlinearity, involving two conden- 
sate fields, the fermion density fluctuations which are 
quadratic in the fermion fields, and the fermion field. 
Since these five fields can mix to produce a sixth, the 
physical process involved is six-wave mixing between the 
boson and fermion fields. Finally, the last term in Eq. (|^), 
r<r = /er</>B(r)[£( )(r,i) + h.c], is a stochastic poten- 
tial whose physical origin are the density fluctuations of 



the vacuum state of the Bose quasi-particles, and has 
the two-point correlation function (r CT (r, t)T a i (r', t')) = 
n B &aa> (r,r',i - t'). 

To proceed we next assume that the fermion beam 
propagates through the condensate at velocity v less 
than the condensate sound velocity c = h^/ATtnsaB /mB- 
Viewing the beam as an impurity traversing the conden- 
sate, the condition v < c implies that it will not create 
incoherent phonon excitations that persist in the con- 
densate after the beam has passed [16). Physically, in 
this limit the fermion beam is accompanied as it prop- 
agates by a virtual phonon cloud that mediates an ef- 
fective instantaneous interaction between the fermions. 
Furthermore, since the phonon excitations are virtual 
the effects of the stochastic potential r CT (r, t) may safely 
be neglected. Thus, we may neglect time retardation 
in the collision term and the stochastic potential in the 
fermionic field Eq. <^j. It then reduces to 
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di) a (r,t) 



(H Fa - nS) + h-4>i4<T'4>o + d 3 r , U <rtr >(r,r f )5h tr >(T f ,t)'ip <r (T f ,t). 



(6) 
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FIG. 1: Schematic diagram for the formation of the two-atom 
bound state in a fermionic beam propagating through a Bose- 
Einstein condensate 



where for the case of a homogeneous condensate 
?7 CTCT '(r, r') is the Yukawa potential & M 



U aa ,{T,T>) 



\ TTk 2 



-V2\T-r'\/l h 



r — r' 



(7) 



and lh = h/^2rnB~gnB is the healing length of the Bose 
condensate. The minus sign in Eq. (Q) indicates that the 
Fermi-Bose coupling induces an attractive effective force 
between fermionic atoms. In the language of nonlinear 
optics the condensate acts as a "nonlinear crystal" for 
the fermions and the induced effective fermion-fermion 
interaction U arT i (r, r', t) plays the role of a spatially non- 
local Kerr nonlinearity for the fermionic field, leading 
to familiar effects such as self-focusing. However, since 
fermionic fields are intrinsically multimode and are not 



amenable to a mean-field description |l7| one cannot ex- 
pect simple, few-mode nonlinear effects such as occur in 
optics or bosonic atom optics, but rather multimode cou- 
pling and many-particle quantum correlations. Here we 
consider the particular case of formation of two-fermion 
bound states, a problem that presents a close analogy to 
the quantum propagation of a two-photon light beam in 
a self- focusing medium jL8j. 

Figure 1 shows a potential nonlinear atom op- 
tics scheme to realize the bound states. We as- 
sume that an ultracold source of fermionic atoms 
generates a beam containing simultaneously two 
atoms. First we consider a spin-polarized beam de- 
scribed by the two-atom quantum state, = 
■j- f d 3 n J d 3 r 2 /(r!,r 2 , 0^(1-1)^2)10}, that then 
propagates at velocity v into a condensate. The func- 
tion /(ri,r 2 ,i) gives the probability amplitude to find 
one atom at ri and the other at r 2 , and is antisymmet- 
ric under atomic exchange. (We ignore the spin label for 
simplicity of notation.) The time evolution of |$(t)) is 
determined by the Schrodinger equation 



ifi 



gjg) 



(8) 



where the total effective Hamiltonian, constructed from 
Eq.©, is 

ffeff = / d 3 rft(r) {H F - (m f ) $(r) 



d A rd A r'U{v, v')^(Y)5n{v')^{v). (9) 



This yields the equation of motion for /(ri, r 2 , t) 



J 



if) 



d/(n,r 2 ,£) 
dt 



= (H F (v 1 ) + H F (v 2 ) + U(r 1 ,r 2 ) + 2 d A r'd A r" [U(r u r') + U(v 2 , r')] |/(r', r", t)\ 2 ) /(n, r 2 , t), 



(10) 



which clearly shows the role of the effective Kerr medium 
on the fermionic field. To determine the condition for 
formation of the two-atom bound states of Eq. jlO| ) 
we transform into the center of mass coordinates R = 
(rj + r 2 )/2 and r = ri — r 2 , and make the ansatz 
/(ri, r 2 , t) = e 'K-R-.Ei/»n/( r ) ] w hich leads finally to the 
eigenvalue problem for the relative motion of two parti- 



cles in a Yukawa potential 



2/x 



W{r) = E r W(r) 



(11) 



where Uo = f 2 mBnB / {kTi 2 ) is the potential strength, 
(i = m F /2 the reduced mass, and E r the energy of rela- 
tive motion. From Eq. dill) , one can see that the spatial 
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width of the energy eigenstates is dictated by the length 
scale l = h 2 /(iiU ). The energy eigenvalues of a Yukawa 
potential have been extensively investigatetd. An accu- 
rate formula has been provided for the number of 
bound states for a state with total angular momentum 
quantum number I, 



(12) 



U=(VZ- y/Zl)Sl + 1, 



where Z = l h /(V2l ),Zi = 0.8399(1 + 2.7359/ + 1.6242Z 2 ) 
and Si = 1.1335(1 + 0.0191/ - 0.001684Z 2 ). Considering 
the requirement of symmetry of the spatial wave func- 
tion for the spin-polarized beam, the angular momentum 
must be I = 1,3,5.... The condition that there exists at 
least one bound state is Z > Z\ . This determines a spa- 
tial range Iq < Ij, = Z^/6.366 for two-atom bound state in 
a Bose condensate, and meanwhile enforces a requirement 
of sufficient potential strength Uq > Ub = 6. 366S 2 /(fJ,lh) 
for binding. 

Possible candidates to observe two-fermion bound 
states include combinations of alkali atoms such as 6 Li- 
7 Li, 6 Li- 23 Na, and 40 K- 39 K. The 40 K- 39 K combination 
appears to be the most promising, due to its very large 
boson-fermion scattering length cibf ~ 52.9 nm and its 
very small boson-boson scattering length ob ~ 0.26 nm. 
In this case, the realization of bound fermionic states re- 
quires a condensate density of tlb ~ 10 13 cm -3 , a value 
achievable with current techniques. 

The basic signature of the formation of a two-fermion 
bound state in the condensate is that it is immune to 
wave packet spreading and can therefore be collected at a 
localized detector placed at the output of the condensate 
as illustrated in Fig. 1. In contrast, if the condensate is 
absent, and the time of flight of the fermions is chosen 
long enough, the initial fermion packet will spread out 
such as to negligibly overlap the detector. The signature 
of the formation of bound states in this nonlinear atom 
optics arrangement is therefore an enhanced detection 
of atoms when the condensate is in place in contrast to 
when it is absent. This remains true even though the 
initial fermion wave packet does not exactly match the 
bound state, as long as it has some reasonable projection 
onto the bound state wave function. 

We remark that our results are readily generalized to 
fermionic atomic beam with two internal spin states. 
In particular, if one initially prepares the beam in a 
spin-singlet state, the two-atom quantum state will have 
the form |$) - ± Jf ^nd^n, r 2 , i)(V>}(rity{(r 2 ) - 

■0j(ri)"(/>j(r2))|O). Then in addition to the attractive 
Yukawa potential the two fermions also experience an 
s-wave interaction, since they are in different spin states. 
Nonetheless two-fermion bound states can still be formed 
under suitable combination of these potentials. In par- 
ticular, the formation of the spatial bound states in the 
spin-singlet case may be useful to create the spatially 



entangled fermionic atomic beam. 

In conclusion, we have derived a nonlinear Schrodinger 
equation for a fermionic beam in a condensate, and 
shown that the condensate acts as an effective nonlinear 
medium for the fermions, leading to formation of a two- 
atom bound state closely analogous to the two-photon 
bound state in self- focusing media |l^]. A clear signa- 
ture of the bound state in a nonlinear atom optics ex- 
periment has also been proposed. Physically, the under- 
standing of fermionic bound states may be important for 
the manipulation of the quantum statistical properties of 
fermionic atomic beams, e.g. changes from antibunched 
and bunched beams, dynamic Cooper pairing, and po- 
tentially the formation of quantum solitons in ultracold 
fermionic atomic beams. 
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